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Abstract 



. In this paper we investigate the encoding of operator quantum error 

correcting codes i.e. subsystem codes. We show that encoding of 
subsystem codes can be reduced to encoding of a related stabilizer 
code making it possible to use all the known results on encoding of 
stabilizer codes. Along the way we also show how Clifford codes can 
q-( be encoded. We also show that gauge qubits can be exploited to reduce 

the encoding complexity 



I> ' Introduction. In this paper we investigate encoding of subsystem codes. 



Our main result is that encoding of a subsystem code can be reduced to the 
q{ , encoding of a related stabilizer code, thereby making use of the previous 

theory on encoding stabilizer codes [2-4]. We shall prove this in two steps. 
\& • First, we shall show that Clifford codes can be encoded using the same 

methods used for stabilizer codes. Secondly, we shall show how these meth- 
yl \ ods can be adapted to encode Clifford subsystem codes. Since subsystem 
codes subsume stabilizer codes, noiseless subsystems and decoherence free 
subspaces, these results imply that we can essentially use the same methods 
to encode all these codes. In fact, while the exact details were not provided, 
it was suggested in [10] that encoding of subsystem codes can be achieved 
by Clifford unitaries. Our treatment is comprehensive and gives proofs for 
all the claims. 

Subsystem codes can potentially lead to simpler error recovery schemes. 
In a similar vein, they can also simplify the encoding process, though per- 
haps not as dramatically 1 . These simplifications have not been investigated 



1 In general, decoding is usually of greater complexity than encoding and for this reason 
it is often neglected in comparison. This parallels the classical case where also the decoding 
is studied much more extensively than encoding. 
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thoroughly, neither have the gains in encoding been fully characterized. Es- 
sentially, these gains are in two forms. In the encoded state there need not 
exist a one to one correspondence between the gauge qubits and the physical 
qubits. However, prior to encoding such a correspondence exists. We can 
exploit this identification between the virtual qubits and the physical qubits 
before encoding to tolerate errors on the gauge qubits, a fact which was 
recognized in [10]. Alternatively, we can optimize the encoding circuits by 
eliminating certain encoding operations. The encoding operations that are 
saved correspond to the encoded operators on the gauge qubits. This is a 
slightly subtle point and will be elaborated at length subsequently. We argue 
that optimizing the encoding circuit for the latter is much more beneficial 
than simply allowing for random initialization of gauge qubits. 

Notation. We shall denote a finite field with q elements by ¥ q . Following 
standard convention we use [[n, k, d]] q for stabilizer codes and [[n, k, r, d]] q 
for subsystem codes. The inner product of two characters of a group N, say 
X and 9, is defined as (x,6)n = J2neN x( n )^( n 1 )- We shall denote 

the center of a group N by Z(N). Given a subgroup N < E,we shall denote 
the centralizer of N in E by Ce(N). Given a matrix A, we consider another 
matrix B obtained from A by column permutation n as being equivalent and 
denote this by B = n A. Often we shall represent the basis of a group by the 
rows of a matrix. In this case we will regard another basis obtained by any 
row operations or permutations as being equivalent and by a slight abuse of 
notation continue to denote B = 7T A. The commutator of two operators A, 
B is defined as [A, B] = AB - BA. 

Encoding Stabilizer Codes [2, 3]. We shall now briefly, review the 
standard form encoding of stabilizer codes, due to Cleve and Gottesman, 
see [2,3]. Recall the Pauli matrix operators 2 , 



X = 



"01" 




" 1 


, Y = 


" -1 " 


, z = 


10 




0-1 


1 



= xz. 



Let V n be the Pauli group on n qubits. An element element e = (— l) c X ai Z bl $ 
•••(g) X a "Z bn in V n , can be mapped to F 2 ,™ by r : V n -> F% n as 

r(e) = (a 1 ,...,a n \b 1 ,...,b n ). (2) 

Given an [[n, k, d\\2 code with stabilizer S, we can associate to S (and 
therefore to the code), a matrix in F^™ k ^ x2n obtained by taking the image 



2 We consider the real version of the Pauli group in this paper. 
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of any set of its generators under the mapping r. We shall refer to this 
matrix as the stabilizer matrix. We shall refer to the stabilizer as well as 
any set of generators as the stabilizer. Additionally, because of the mapping 
r, we shall refer to the stabilizer matrix or any matrix obtained from it by 
row reduction or column permutations also as the stabilizer. The stabilizer 
matrix can be put in the so-called "standard form", see [2,3]. This form 
also allows us to compute the encoded operators for the stabilizer code. 
Recall that the encoded operators allow us to perform computations on the 
encoded data without having to decode the data and then compute. 

Definition 1 (Encoded operators). Given a [[n,k,d]]2 stabilizer code with 
stabilizer S, let Xi, Z{ for 1 < i < k be a set of 2k linearly independent 
operators in C-p n {S) \ SZ(V n ). The operators Xi, Zi are said to be encoded 
operators for the code if they satisfy the following requirements, 
i) {XuXj} = 
it) [Z i ,Zj] = Q __ 
Hi) [Xi, Zj] = 25ijXiZi 

The operators Xi and Zj are referred to as encoded or logical X and Z 
operators on the ith and j'th logical qubits, respectively. The choice of which 
of the 2k linearly independent elements of C-p n (S)\SZ(V n ) we choose to call 
encoded X operators and Z operators is arbitrary; as long as the generators 
satisfy the conditions above, any choice is valid. Different choices lead to 
different sets of encoded logical states; alternatively, a different orthonormal 
basis for the codespace. 

Lemma 2 (Standard form of stabilizer matrix [2,3]). Up to a permutation 
TT, the stabilizer matrix of an [[n, k, d]]2 code can be put in the following form, 







A 1 




A 2 
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D 
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C 
E 



while the associated encoded operators can be derived as 



z ' 




"00 


A\ I k ~ 
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7T 


[ E 1 I k 


C l _ 



(3) 



(4) 



Remark 3. Encoding using essentially same ideas is possible even if the 
identity matrices (I s > in the stabilizer matrix or in the encoded operators) 
are replaced by upper triangular matrices. 

The standard form of the stabilizer matrix prompts us to distinguish 
between two types of the generators for the stabilizer as they affect the en- 
coding in different ways (although it can be shown that they are of equivalent 
complexity). 
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Definition 4 (Primary generators). A generator Gi = (a±, . . . , a n \bi, . . . , b n ) 
with at least one nonzero ai is called a primary generator. 

In other words, primary generators contain at least one X 01 Y operator 
on some qubit. The primary generators determine to a large extent the 
complexity of the encoding circuit along with the encoded X operators. 
The operators X are also called seed generators and they also figure in the 
encoding circuit. The encoded Z operators do not. 

Definition 5 (Secondary generators). A generator of the form (0, . . . , 0|&i, . . . , b n ) 
is called secondary generator. 

In the standard form encoding, the complexity of the encoded X oper- 
ators is determined by the secondary generators. Therefore they indirectly 
contribute 3 to the complexity of encoding. 

We mentioned earlier that different choices of the encoded operators 
amounts to choosing different orthonormal basis for the codespace. How- 
ever, the choice in Lemma 2 is particularly suitable for encoding. We can 
represent our input in the form |0)® \a± . . . atk) which allows us to make 
the identification that |0}® is mapped to |o), the logical all zero code word. 
This state is precisely the state stabilized by the stabilizer generators and 
logical Z operators, (which in Lemma 2 can be seen to be consisting of only 
Z operators). Given the stabilizer matrix in the standard form and the 
encoded operators as in Lemma 2, the encoding circuit is given as follows. 

Lemma 6 (Standard form encoding of stabilizer codes [2,3]). Let S be the 

stabilizer matrix of an [[n, k, d\\2 stabilizer code in the standard form i.e., 
as in equation (3). Let Gi denote the ith primary generator of S and Xj 
denote the jth encoded X operator as in equation (4). Then these operators 
are in the form 4 

Gi = (0,0, . . . , l,a i+ i, . . . ,a n \bi, . . . ,b s ',0, . . . ,0,b n - k+ i, ■ ■ ■ ,b n ), 
Xj = (0, . . . ,0,cy+i, . . . ,c„_ fc 0, ... ,0, 1 = c n - k+j ,0, . . . ,0\di, . . .,d s ',0, . . . ,0). 

To encode the stabilizer code we implement the following circuits correspond- 
ing to each of the primary generators and the encoded operators. The gen- 
erator Gi is implemented after Gj+i. The encoded operators precede the 

3 Indirect because the submatrix E, figures in both the secondary generators, see equa- 
tion (3), and also the encoded X operators, see equation (4). 

4 We allow some freedom in the primary generators, in that instead of I a > in equation (3), 
we allow it be an upper triangular matrix also. 
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primary generators in their implementation but we can implement Xj be- 
fore or after Xj +i . 



\0} 1 
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X a ™-k+l Z bn ' k + 1 
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To encode a stabilizer code, we first put the stabilizer matrix in the 
standard form, then implement the seed generators i.e., the encoded X 
operators, followed by the primary generators i = s' to i = 1 as per Lemma 6. 
The complexity of encoding the ith primary generator is at most n — i 
two qubit gates and one H gate. The complexity of encoding an encoded 
operator is at most n — k — s' CNOT gates. This means the complexity of 
standard form encoding is upper bounded by (2n — 1 — k — s')s'/2 two qubit 
gates and s' Hadamard gates; 0(n(n — k)) gates. Perhaps an example will 
help at this juncture. 
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Example 7. Let us consider the [[5,1,3]] code with following stabilizer. 



S = 



X I X X X 

I X Z X Y 

Z I Z Z Z 

I Z Y Z X 
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The associated stabilizer matrix is given by 



S = 



Writing S in standard form we get 



S = 



The encoded operators for this code are 
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The stabilizer matrix has three primary generators. By Lemma 6 the encod- 
ing circuit is given by 
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|0> 
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|o> —\x 



-t- 



7=1- 
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Z 



X G 3 G 2 d 

Figure 1: Encoding for the [[5, 1,3]] code 
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Remark 8. At this point we want to highlight that in Lemma 6, we let 
primary generators to be in the upper triangular form instead of the stan- 
dard form given in Lemma 2, because of which the primary generators were 
required to be implemented in a particular order. If however, we had them 
strictly in the standard form of Lemma 2 then any order is possible. For 
instance, implementing the generators in the reverse order for the [[5,1,3]] 
will give the following circuit. However, note that some additional Z gates 



10) 
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X G\ G2 G3 

Figure 2: Alternative encoding for the [[5, 1,3]] code 

are present in this implementation for G\, while some of the Z gates asso- 
ciated to G\ and G2 are redundant. The reduced circuit will be given as: 
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Figure 3: Alternative encoding for the [[5, 1, 3]] code with redundant Z gates 
removed 

Other orderings of the primary generators are possible, but we must be 
careful to account for the Z gates that do not act on the |0) state directly. 



Encoding Clifford Codes. In this section, we show that a Clifford code 
can be encoded using its stabilizer and therefore the methods used for en- 
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coding stabilizer codes are applicable. We briefly recapitulate some facts 
about Clifford subsystem codes. Let E be an abstract error group i.e., it is 
a finite group with a faithful irreducible unitary representation p of degree 
\E : Z(E)\ 1 / 2 . Denote by <f>, the irreducible character afforded by p. Let ./V 
be a normal subgroup of E. Further, let x be- an irreducible character x of 
N such that ((()n,x)n > 0. Then the Clifford code defined by (E, p, N, x) is 
the image of the orthogonal projector 



Under certain conditions we can construct a subsystem code from the 
Clifford code, in particular when E is the extraspecial p-group, the Clifford 
code C has a tensor product decomposition 5 as C = A <g> B, where B is an 
irreducible CiV-module, A is an irreducible CL- module and L = Ce(N). In 
this case we can encode information only into the subsystem A, while the 
co-subsystem B provides additional protection. When encoded this way we 
say C is a Clifford subsystem code. The normal subgroup N consists of 
all errors in E that act trivially on A. It is also called the gauge group of 
the subsystem code. Our main goal will be to show how to encode into the 
subsystem A. Therefore, our interest will center on the projectors for the 
Clifford code and the subsystem code and not so much on the parameters 
of the codes themselves. 

An alternate projector for a Clifford code with data (E, p, N, x) can be 
defined in terms of Z(N), the center of N. The proof of this can be found 
in [6, Theorem 6]. This projector is given as 



where ip is an irreducible character of Z(N), that satisfies (x { Z(N))(x) = 
x(l)<p(x). In this case Q can be thought of as a stabilizer code in the sense 
of [1] i.e. 



In addition to the assumption that the error group is an extraspecial p-group 
we also assume that Z(E) < N. The inclusion of the center of E does not 
change the code but helps in analysis. Thus we have the following lemma. 

5 Strictly speaking the equality should be replaced by an isomorphism. 




(5) 




(6) 



pirn) = (p(m) \ip) for any m in Z(N). 



(7) 
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Lemma 9. Let (E, p, N, x) be the data of a Clifford code and ip an irreducible 
character of Z(N), the center of N, satisfying (x i Z(N))(x) = x(l)y(x). 
If E is an extraspecial p-group, then for all n in Z(N), ip(n) G {C k | C = 
e j2nk/ P ^Q < k < p}. Further, if Z(E) < N, then for any n G Z(N), we have 
v(n- l )p{n)ep{Z{N)). 

Proof. First we note that the irreducibilty of p implies that for any z in Z{E) 
we have p(z) = loI for some u G C by Schur's lemma. The assumption that 
E is an extraspecial p-group forces to G {£ fc | < k < p} where £ = e ?27r//p . 
This is because ^(E 1 )! = p for extraspecial p-groups. Secondly, we observe 
that ip is an irreducible additive character of Z(N) (an abelian subgroup 
of an extraspecial p-group) which implies that we must have <p{ri) = C} for 
some < I < p, [8]. Together these observations imply that we can assume 
^(n- 1 )/ = ( l I = p(z) for some < I < p and z G Z(E). Since Z(E) < N, 
it follows that Z(E) < Z(N) and v( n - l )p{n) is in p(Z(N)). □ 

Our goal is to use the stabilizer of Q for encoding and as a first step 
we will show that it can be computed from Z(N). The usefulness of such a 
projector is that it obviates the need to know the character ip. Let S < p(E) 
be the stabilizer of Q. Then we claim that S is given as 

S = {(^(n _1 )p(n) | n G Z(N)}. 

We claim that S can be used for encoding the associated Clifford code. Then 
we will show how the encoding circuit of the Clifford code is to be modified 
so that we can encode the subsystem code derived from the Clifford code. 

Theorem 10. Let Q be a Clifford code with the data (E, p, N, \) and ip a 
constituent of the restriction of x to Z = Z(N). Let E be an extraspecial 
p-group and Z{E) < N and 

S={<p(n- 1 )p(n)\n£Z(N)} and P = ^ £ s. (8) 

Then S is the stabilizer of Q and Im P = Q. 
Proof. We will show this in a series of steps. 

1) First we will show that S < p{Z). By Lemma 9 we know that (p(n~ 1 )p(n) 
is in p(Z), therefore S C p(Z). For any two elements n±,U2 G Z, we 
have si = <p{n~l l ) p{n\) , S2 = (^(n^ 1 )p(ri2) G S and we can verify that 
s^ 1 s 2 = l f(ni)p(n^ 1 )p(n 2 1 )p(n 2 ) = •p{n' 2 1 ni)p(n^ 1 n2) G S, as p(n^ 1 n 2 ) 
is in p(Z). Hence S < p(Z). 
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2) Now we show that S fixes Q. Let s G S and |^} G Q. Then s = 
Lp(n~ l )p(n) for some n G The action of s on is given as s \ip) = 
ip(n~ 1 )p(n) \ip) = ip(n~ 1 )(p(n) \tp) = \tp), in other words S fixes Q. 

3) Next, we show that \S\ = \Z\/\Z(E)\. If two elements rt\ and n<i in Z 
map to the same element in S, then (^(n ] " 1 )p(ni) = ipiri^ 1 ) p(ri2) ■, that is 
p(^2) = t p(n± 1 n2)p(ni). From Lemma 9 it follows that p{n 2 ) = CV( n i) 
for some < I < p. Since p(Z(E)) = { e j27rk / p I | < k < p}, we must 
have n 2 = zm for some z G Z(£7). Thus, \S\ = \Z\/\Z(E)\. 

4) Let T be a traversal of Z(E) in Z, then every element in Z can be written 
as zt for some z G Z(.E) and i G T. From step 3) we can see that all 
elements in a coset of Z(E) in Z map to the same element in S, therefore, 

S = {<p(t- 1 ) P (t)\teT}. 

Recall that a projector for Q is given by 

' ' neZ 

But we know from step 3) that if z G Z(E), then <^(n -1 )/9(n) = v?((2;n)~ 1 )/9(zn). 
So we can simplify P' as 

p ' = E ^m*), 

1 1 teT zeZ(E) 



i^i ^ 



teT 
s = P. 



ses 

Thus the projector defined by S 1 is precisely the same as P' and P is also 
a projector for Q. 

From step 3) it is clear that S n Z(E) = {1} and by [5, Lemma 10], S is 
a closed subgroup of E. By [5, Lemma 9], Im P = Q is a stabilizer code. 
Hence 5 is the stabilizer of Q. □ 
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Corollary 11. Let Q be an [[n, k,r, d]] Clifford subsystem code and S its 
stabilizer. Let 



Then P is a projector for the subsystem code i.e. Q = Im P. 

Proof. By [7, Theorem 4], we know that an [[n, k, r, d]] Clifford subsystem 
code is derived from a Clifford code with data (E, p, N, This construction 
assumes that E is an extraspecial p-group and Z(E) < N < E. Since as 
subspaces the Clifford code and subsystem code are identical, by Theorem 10 
we conclude that the projector defined from the stabilizer of the subspace is 
also a projector for the subsystem code. □ 

Theorem 10 shows that any Clifford code can be encoded using its sta- 
bilizer. As to a subsystem code, while Corollary 11 shows that there exists 
a projector that can be defined from its stabilizer, it is not clear how to 
use it so that one respects the subsystem structure during encoding. More 
precisely, how do we use the projector defined in Corollary 11 to encode into 
the information carrying subsystem A and not the gauge subsystem. This 
will be the focus of the next section. 

Encoding Subsystem Codes. For ease of presentation and clarity hence- 
forth we will focus on binary codes, though the results can be extended to 
nonbinary alphabet using methods similar to stabilizer codes, see [4]. The- 
orem 10 shows that in order to encode Clifford codes we can use a projector 
derived from the underlying stabilizer to project onto the codespace. But 
in case of Clifford subsystem codes we know that Q = A® B and the infor- 
mation is to be actually encoded in A. Hence, it is not sufficient to merely 
project onto Q, we must also show that we encode into A when we encode 
using the projector defined in Corollary 11. 

Let us clarify what we mean by encoding the information in A and not 
in B. Suppose that P maps |0) to A <8> |0) B and |1) to A <8> |1) B - Then 
the information is actually encoded into B. Since the gauge group acts 
nontrivially on B, this particular encoding does not protect information. Of 
course a subsystem code should not encode (only) into B, but we have to 
show that the projector defined by P s does not do that. 

We need the following result on the structure of the gauge group and the 
encoded operators of a subsystem code. Poulin [9] proved a useful result on 
the structure of the gauge group and the encoded operators of the subsystem 




(9) 
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code. But first a little notation. A basis for V n is Xi, Zi, 1 < i < n, where 
Xi and Zi are given as 

n n 

X i = (g)X^ and Z i = (g)Z*y. 
j"=i j"=i 

They satisfy the relations [X u Xj\ = = [Z^ZA; [X u Zj] = 25 ij X i Z j . 
However, we can choose other generating sets {xj, zi | 1 < i < n} for V n that 
satisfy similar commutation relations i.e., [xi,Xj] = = [zi,Zj] and [xj,Zj] = 
25{jXiZj. These operators may act nontrivially on many qubits. We often 
refer to the pair of operators Xi,Zi that satisfy the commutation relations 
similar to the Pauli operators as a hyperbolic pair. Given an [[n, k, r, d]] 
code we could view the state space of the physical n qubits as that of n 
virtual qubits on which these Xi, Zi act as X and Z operators. In particular 
k of these virtual qubits are the logical qubits and r of them gauge qubits. 
The usefulness of these operators is that we can specify the structure of the 
stabilizer, the gauge group and the encoded operators. The following lemma 
makes this specification precise. 

Lemma 12. Let Q be an [[n, k, r, c£]] 2 subsystem code with gauge group, G 
and stabilizer S. Denote the encoded operators by Xi, Zi, 1 < i < k, where 
[Xi,Xj] = = [Zi, Zj];[Xi, Zj] = 25ijXiZj. Then there exist operators 
{xi, Zi G V n | 1 < i < n} such that 

i) S = {zi,z 2 , ■ ■ -,z s ), 

ii) G = (S, z s+ i,x s+ i, . . . , z s+r , x 3 +r, Z{V n )), 
iH) Gp n {S) = (G^uZu ...... .,X k ,Z k ), 

iv) Xi = x s+r+ i and Zi = z a +r+i> 1 — * — k, 

where [zi,Zj] = [xi,Xj] = 0; [xj,Zj] = 25ijXiZi. Further, S defines an [[n,k + 
r]] stabilizer code encoding into the same space as the subsystem code and its 
encoded operators are given by {x s+ i, z s+ \, . . . , x s+r , z s+r , X\, Z\,..., X k , Z^} 

Proof. See [9] for proof on the structure of the groups. Let Q = A® B, then 
dim A = 2 k and dimB = 2 r . From Corollary 11 we know that the projector 
defined by S also projects onto Q (which is 2 fc+r -dimensional) and therefore 
it defines an [[n, k + r]] stabilizer code. From the definition of the operators 
Xi, Zi and Xi, Zi and the fact that 

C-Pn{S) = (S, x s+ i, z s+ i, . . . , x s+r , z s+r Xi, Z\, . . . , Xk, Z k , Z{V n )) 

we see that x« , Z{ , for s + 1 < i < r act like encoded operators on the gauge 
qubits, while Xi, Zi continue to be the encoded operators on the information 
qubits. Together they exhaust the set of 2(k + r) encoded operators of the 
[[n, k + r]] stabilizer code. □ 
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We observe that the logical operators of the subsystem code are also 
logical operators for the underlying stabilizer code. So if the stabilizer code 
and the subsystem code have the same logical all zero state, then Lemma 12 
suggests that in order to encode the subsystem code, we can treat it as 
stabilizer code and use the same techniques to encode. If the logical all zero 
code word was the same for both the codes, then because they have the same 
logical operators we can encode any given input to the same logical state 
in both cases. Using linearity we could then encode any arbitrary state. 
Encoding the all zero state seems to be the key. Now, even in the case of 
the stabilizer codes, there is no unique all zero logical state. There are many 
possible choices. The reader can refer to the appendix for examples. Given 
the encoded operators it is easy to define the logical all zero state as the 
following definition shows: 

Definition 13. A logical all zero state of an [[n, k, r, d]] subsystem code is 
any state that is fixed by its stabilizer and k logical Z operators. 

This definition is valid in case of stabilizer codes also. This definition 
might appear a little circular. After all, we seem to have assumed the 
definition of the logical Z operators. Actually, this is a legitimate definition 
because, depending on the choice of our logical operators, we can have many 
choices of the logical all zero state. In case of the subsystem codes, this 
definition implies that the logical all zero state is fixed by n — r operators, 
consequently it can be any state in that 2 r -dimensional subspace. If we 
consider the [[n, k + r]] stabilizer code that is associated to the subsystem 
code, then its logical zero is additionally fixed by r more operators. So 
any logical zero of the stabilizer code is also a logical all zero state of the 
subsystem code. It follows that if we know how to encode the stabilizer 
code's logical all zero, we know how to encode the subsystem code. We are 
interested in more than merely encoding the subsystem code of course. We 
also want to leverage the gauge qubits to simplify and/or make the encoding 
process more robust. Perhaps a few examples will clarify the ideas. 

Illustrative Examples. Consider the following [[4, 1, 1,2]]2 subsystem 
code, with the gauge group G, stabilizer S and encoded operators given 
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The encoded operators of this code are given by 
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The associated [[4, 2]] stabilizer code has the following encoded operators. 
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It will be observed that the encoded X operators of [[4, 2]] are in a form 
convenient for encoding. We treat the [[4,1,1,2]] code as [[4,2]] code and 
encode it as in Figure 4. The gauge qubits are permitted to be in any state. 



10} 
l<7> 

|0> 

Figure 4: Encoding the [[4, 1, 1, 2]] code (Gauge qubits can be in any state) 
Assuming g = a |0} + 6 11), the logical states up to a normalizing constant 
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are 

|0> = o(|0000) + |llll))+6(|0101) + |1010», 
|T) = a(|0011) + |1100» + &(|0110) + |1001». 

It can be easily verified that S stabilizes the above state and while the gauge 
group acts in a nontrivial fashion, the resulting states are still orthogonal. 
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In this example we have encoded as if we were encoding the [[4, 2]] code. 
Prior to encoding the gauge qubits can be identified with physical qubits. 
After the encoding however such a correspondence between the physical 
qubits and gauge qubits does not necessarily exist in a nontrivial subsystem 
code. Since the encoded operators of the subsystem code are also encoded 
operators for the stabilizer code, we are guaranteed that the information is 
not encoded into the gauge subsystem. 

As the state of gauge qubits is of no consequence, we can initialize them 
to any state. Alternatively, if we initialized them to zero, we can simplify 
the circuit as shown in Figure 5. 
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Figure 5: Encoding the [[4, 1, 1, 2]] code (Gauge qubits initialized to zero) 

The encoded states in this case are (again, the normalization factors are 
ignored) 

|0> = |0000) + |1111} , 
|T) = |0011) + |1100} . 

The benefit with respect to the previous version is that at the cost of ini- 
tializing the gauge qubits, we have been able to get rid of all the encoded 
operators associated with them. This seems to be a better option than ran- 
domly initializing the gauge qubits. Because it is certainly easier to prepare 
them in a known state like 1 0) , rather than implement a series of controlled 
gates depending on the encoded operators associated with those qubits. 

At this point we might ask if it is possible to get both the benefits of 
random initialization of the gauge qubits as well as avoid implementing the 
encoded operators associated with them. To answer this question let us look 
a little more closely at the previous two encoding circuits for the subsystem 
codes. We can see from them that it will not work in general. Let us see 
why. If we initialize the gauge qubit to |1) instead of |0) in the encoding 
given in Figure 5, then the encoded state is 

|0) = |0100) + |1011} , 
|T> = |0111) + |1000} . 
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Both these states are not stabilized by S, indicating that these states are 
not in the code space. 

In general, an encoding circuit where it is simultaneously possible initial- 
ize the gauge qubits to random states and also avoid the encoded operators 
is likely to be having more complex primary generators. For instance, let us 
consider the following [[4, 1, 1, 2]] subsystem code: 
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The encoded operators of this code are given by 



L = 
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The associated [[4, 2]] stabilizer code has the following encoded operators. 
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The encoding circuit for this code is given by 
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Figure 6: Encoding [[4, 1, 1, 2]] code (Encoded operators for the gauge qubits 
are trivial and gauge qubits can be initialized to random states) 



In this particular case, the gauge qubits (as well as the information 
qubits) do not require any additional encoding circuitry. In this case we 
can initialize the gauge qubits to any state we want. But, the reader would 
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have observed we did not altogether end up with a simpler circuit. The 
primary generators are two as against one and the complexity of the encoded 
operators has been shifted to them. So even though we were able to get 
rid of the encoded operator on the gauge qubit and also get the benefit 
of initializing it to a random state, this is still more complex compared to 
either of encoders in Figures 4 and 5. Our contention is that it is better 
to initialize the gague qubits to zero state and not implement the encoded 
operators associated to them. 

Encoding Subsystem Codes by Standard Form Method. The pre- 
vious two examples might lead us to conclude that we can take the stabilizer 
of the given subsystem code and form the encoded operators by reducing 
the stablizer to its standard form and encode as if it were a stabilizer code. 
However, there are certain subtle points to be kept in mind. When we form 
the encoded operators we get k + r encoded operators; we cannot from the 
stabilizer alone conclude which are the encoded operators on the informa- 
tion qubits and which on the gauge qubits. Put differently, these operators 
belong to the space C-p n (S) \ S = GC-p n (G) \ SZ(V n ). It is not guaranteed 
that they are entirely in Cp n (G) i.e., we cannot say if they act as encoded 
operators on the logical qubits. This implies that in general all these oper- 
ators act nontrivially on both A and B. Consequently, we must be careful 
in choosing the encoded operators and the gauge group must be taken into 
account. We give two slightly different methods for encoding subsystem 
codes. The difference between the two methods is subtle. Both methods 
require the gauge qubits to be initialized to zero. In the second method 
(see Algorithm 2) however, we can avoid the encoded operators associated 
to them. Under certain circumstances, we can also permit initialization to 
random states. 

Correctness of Algorithm 1. Since stabilizer Sa > S, the space 
stabilized by Sa is a subspace of the A <g) B, the subspace stabilized by 
S. As ISaI/I^I = 2 r , the dimension of the subspace stabilized by Sa is 
2 k+r /2 r = 2 k . Additionally, the generators z s +i, ■ ■ ■ , z s + r act trivially on A. 
The encoded operators as computed in the algorithm act nontrivially on A 
and give 2 k orthogonal states; thus we are assured that the information is 
encoded into A. 

Let us encode the [[9,1,4,3]] Bacon-Shor code using the method just 
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Algorithm 1 ENCODING SUBSYSTEM CODES - STANDARD FORM METHOD 1 

Require: Gauge group, G = (S,x s +i, z s +i, ■ ■ ■ ,x s+r , z s+r ,H) and stabi- 
lizer, S = (z±, . . . , z n -k-r) °f the [[n, k, r, d\] subsystem code. 
Ensure: [x u Xj] = [z^Zj] = 0; [x u Zj] = 2x i z i 5 ij 

1: Form Sa = (S, z s+ ±, . . . , z s+r ), where s = n — k — r 
2: Compute the standard form of Sa as per Lemma 2 
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3: Compute the encoded operators X\, . . . , as 
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4: Encode using the primary generators of Sa and Xi as encoded operators, 
see Lemma 6; all the other (n — k) qubits are initialized to |0). 



proposed. The stabilizer and the gauge group are given 6 by 
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We do not show the identity. 
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Let us form Sa by augmenting S with G z . Then 
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The encoded X and Z operators are XfXgXg and Z\Z±Zt, respectively. 
After putting Sa in the standard form, and encoder for this code is given in 
Figure 7. 
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Figure 7: Encoder for the [[9,1,4,3]] code. This is also an encoder for the 
[[9,1,3]] code. 

If on the other hand we had formed Sa by adding G x instead, then Sa 
would have been 
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The encoded operators remain the same. In this case the encoding circuit 
is given in Figure 8. 
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Figure 8: Encoder for the [[9, 1,4,3]] code with fewer CNOT gates. 

The circuit in Figure 8 has fewer CNOT gates, though the number of 
single qubit gates has increased. Since we expect the implementation of the 
CNOT gate to be more complex than the H gate, this might be a better 
choice. In any case, this demonstrates that by exploiting the gauge qubits 
one can find ways to reduce the complexity of encoding circuit. 

The gauge qubits provide a great degree of freedom in encoding. We 
consider the following variant on standard form encoding, where we try to 
minimize the the number of primary generators. This is not guaranteed to 
reduce the overall complexity, since that is determined by both the primary 
generators and the encoded operators. Fewer primary generators might 
usually imply encoded operators with larger complexity. In fact we have 
already seen, that in the case of [[9, 1,4, 3]]2 code that a larger number of 
primary generators does not necessarily imply higher complexity. However, 
it has the potential for lower complexity. 
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Algorithm 2 ENCODING SUBSYSTEM CODES - STANDARD FORM METHOD 2 

Require: Gauge group, G = (S,x s +i, z s +i, ■ ■ ■ ,x s+r , z s+r ,H) and stabi- 
lizer, S = (z±, . . . , z n -k-r) °f the [[n, k, r, d\] subsystem code. 
Ensure: [x u Xj] = [z^Zj] = 0; [x u Zj] = 2x i z i 5 ij 



1: Compute the standard form of S as per Lemma 2 







A 1 A 2 
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E 



2: Form Sa = (S, z s+ ±, . . . , z s+r ), where s = n — k — r 
3: Compute the standard form of Sa as per Lemma 2 
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4: Compute the encoded operators X\, . . . , as 
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5: Encode using the primary generators of S and Xi as encoded operators, 
accounting for tt\ and 7T2, see Lemma 6; all the other (n — /c) qubits are 
initialized to |0). 



The main difference in the second method comes in lines 1 and 5. We 
encode using the primary generators of the stabilizer of the subsystem code 
instead of the augmented stabilizer. The encoded operators however remain 
the same as before. 

Correctness of Algorithm 2. The correctness of this method lies in 
the observation we made earlier (see discussion following Definition 13), that 
any logical all zero state of the stabilizer code is also a logical all zero of the 
subsystem code and the fact that both share the encoded operators on the 
encoded qubits. 

Remark 14. The permutation iT2 in Algorithm 2 can be restricted to the 
last n — s' columns, since while adjoining the additional r generators to S, 
we could take it to be in the standard form. 

The encoded operators are given modulo the elements of the gauge group 
as in Algorithm 1, which implies that the their action might be nontrivial 
on the gauge qubits. The benefit of the second method is when S and Sa 
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have different number of primary generators. The following aspects of both 
the methods are worth highlighting. 

1) The gauge qubits must be initialized to |0) in both methods. 

2) In Algorithm 1, the number of primary generators of S and Sa can 
be different leading to a potential increase in complexity compared to 
encoding with S. 

3) In both methods, the encoded operators as computed are modulo Sa- 
Consequently, the encoded operators might act nontrivially on the gauge 
qubits. 

Encoding Subsystem Codes by Conjugation Method. The other 
benefit of subsystem codes is the random initialization of the gauge qubits. 
We now give circuits where we can encode the subsystem codes to realize 
this benefit. But instead of using the standard form method we will use 
the conjugation method proposed by Grassl et a/., [4] for stabilizer codes. 
After briefly reviewing this method we shall show how it can be modified 
for encoding subsystem codes. 

The conjugation encoding method can be understood as follows. It is 
based on the idea that the Clifford group acts transitively on the Pauli 
error group. It is possible to transform the stabilizer matrix of any [[n, k, d]] 
stabilizer code into the matrix (00|/ n _fcO). For a code with this stabilizer 

matrix the encoding is trivial. We simply map to jO}® \tp) . The 
associated encoded X and Z operators are given by (OI^OO) and (00|07fe) 
respectively. Here we give a sketch of the method for the binary case, the 
reader can refer to [4] for details. Assume that the stabilizer matrix is 
given by S. Then we shall transform it into (00|/ n _fc0) using the following 
sequence of operations. 

(X\Z) i ► (/ n _ fc 0|0) -> (00|J n _ fc 0). (10) 

This can be accomplished through the action of H = [ \ 2i ] , P = [J?] an d 
CNOT gates on the Pauli group under conjugation. The H gate acting on 
the ith qubit on (a±, . . . ,a n \bi, . . . , b n ) transforms it as 

H ■ 

(ai,...,a n \bi,...,b n ) ^4 (oi, . . . , b ; , . . . , a n \h, . . . , a i; . . . , b n ). (11) 

These modified entries have been highlighted for convenience. The phase 
gate P on the ith qubit transforms (ai, . . . , a n \bi, . . . , b n ) as 

p. 

(oi, . . . , a n |6i, . . . , b n ) i-4 (ai, . . . , a ; , . . . , a n |6i, . . . , a ; + b ; , . . . , b n ). (12) 
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We denote the CNOT gate with the control on the ith qubit and the tar- 
get on the jth qubit by CNOTV. The action of the CNOT^ gate on 
(ai, . . . , a n \bi, . . . , b n ) is to transform it to 

(ai, . . . ,aj_i,aj + ai,o J+ i . . . ,o n |6i, . . . , 6»_i, bi + bj, b i+ i, . . . ,a n ). (13) 

Note that the jth entry is changed in the X part while the ith entry is 
changed in the Z part. For example, consider 

(1,0,0,1,0|0,1,1,0,0) ^ (1,0,0,0,0(0,1,1,0,0), 
(1,0,0,1,0|0,1,1,1,0) (1,0,0,0,0|1,1,1,1,0). 

Based on the action of these three gates we have the following lemmas to 
transform error operators. 

Lemma 15. Assume that we have a error operator of the form (a±, . . . , a n \b\, . 
Then we apply the following gates on the ith qubit to transform the stabilizer, 
transforming (aj,6j) to (a, (3) as per the following table. 





Gate 


(a, 13) 


(0,0) 


I 


(0,0) 


(0,1) 


H 


(1,0) 


(1,0) 


I 


(1,0) 


(1,1) 


P 


(1,0) 



Let x denote 1 + x mod 2, then the transformation to (a\ , . . . , a n |0, . . . , 0) is 
achieved by 

n 
i=l 

For example, consider the following generator (1,0, 0, 1,0|0, 1, 1, 1,0). 
This can be transformed to (1, 1, 1, 1, 0|0, 0, 0, 0, 0) by the application of 
I®H®H®P®I. 

Lemma 16. Let e be an error operator of the form (pi, . . . , ctj = 1, . . . , a n |0, . . . 
Then e can be transformed to (0, . . . , 0, a« = 1, 0, . . . , 1 0, . . . , 0) by 

n 

Yl [CNOT^'p. 
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As an example consider (1, 1, 1, 1, 0|0, 0, 0, 0, 0), this can be transformed 
to (0,1,0,0,0(0,0,0,0,0) by 

CNOT 2 ' 1 • CNOT 2 ' 3 • CNOT 2 ' 4 . 

The first step involves making the Z portion of the stabilizer matrix all 
zeros. This is achieved by single qubit operations consisting of H and P 
performed on each row one by one. 

Note that we must also modify the other rows of the stabilizer matrix 
according to the action of the gates applied. 

Once we have a row of stabilizer matrix in the form (a|0), where a is 
nonozero we can transform it to the form (0, ...,0, Oj = 1,0, ...,0|0) by 
using CNOT gates. Thus it is easy to transform (X\Z) to (/ n -fe0|0) using 
CNOT, P and H gates. The final transformation to (0|/ n _fc0) is achieved by 
using H gates on the first n — k qubits. At this point the stabilizer matrix 
has been transformed to a trivial stabilizer matrix which stabilizes the state 
|0}® n * • The encoded operators are (0J fe |0) and (0|0/ fc ). Let T be the 
sequence of gates applied to transform the stabilizer matrix to the trivial 
stabilizer matrix. Then T applied in the reverse order to jO}® gives 
the encoding circuit for the stabilizer code. 

Now we shall use the conjugation method to encode the subsystem codes. 
The main difference is that instead of considering just the stabilizer we need 
to consider the entire gauge group. Let the gauge group be G = (S, Gz, Gx), 
where Gz = (z s +i, • • • , z s+r ), and Gx = {x s +i, ■ ■ ■ , %s+r)- The idea is to 
transform the gauge group as follows. 
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(14) 



At this point the gauge group has been transformed to a group with trivial 
stabilizer and trivial encoded operators for the gauge qubits and the encoded 
qubits. The sequence of gates required to achieve this transformation in the 
reverse order will encode the state |0)® \<p) The state \4>) corresponds to 
the gauge qubits and it can be initialized to any state, while corresponds 
to the input. 
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Algorithm 3 ENCODING SUBSYSTEM CODES - CONJUGATION METHOD 

Require: Gauge group, G = (S,Gz,Gx), where Gz = (z 8 +i, 



and Gx = (x s +i, ■ ■ ■ ,x s + r ) and stabilizer, S 
[[n, k, r, d]] subsystem code. 
Ensure: [xi,xj] = [zi,Zj] = 0; [x^Zj] 

1: Assume that G is the following form 



(z 1 , . . .,z n _ k _ r ) of the 



2>XiZ{ S{j 



G 



S 



G z 



G x 



2 
3 
4 
5 
6 
7 
8 
9 

10: 
11 
12 



for all i = 1 to s + r do 

Transform z-i to z\ = (a±, . . . , a n \0, . . . , 0) using Lemma 15 
Transform z[ to (0, . . . , dj = 1, . . . , 1 0) using Lemma 16 
Perform Gaussian elimination on column i for rows j > i 
end for 

Apply H gate on each qubit i = 1 to i = s + r 

for alH = s + 1 to s + r do 

Transform X{ to x\ = (ai, . . . , a n \0, . . . , 0) using Lemma 15 
Transform x\ to (0, . . . , ctj = 1, . . . , 1 0) using Lemma 16 
Perform Gaussian elimination on column i for rows j > i 

end for 



In the above algorithm, we assume that whenever a row is transformed 
according to Lemma 15 or 16, all the other rows are also transformed ac- 
cording to the transformation applied. 

Correctness of Algorithm 3. The correctness of the algorithm is 
straightforward. As G has full rank of n — k + r, for each row of G, we will 
be able to find some nonzero pair (a, b) so that the the transformation in 
lines 2-6 can be achieved. When S and Gz are in the form (0|/ s+r 0), the 
rows in Gx are in the form 

[0 A B\0 D ]. 

The zero columns of Gx are consequence of the requirement to satsify the 
commutation relations with (transformed) S and Gz- For instance, The 
first n — k — r are all zero because they must commute with (0|/ s 0), the 
elements of the transformed stabilizer. The submatrix A must have rank r, 
otherwise at this point one of the rows of Gx commutes with all the rows of 
Gz and the condition that we have there are r hyperbolic pairs is violated. 
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It is possible therefore to transform A to the form (0/ r 0|0). It cannot be 
any other form because then we would not have the r hyperbolic pairs. The 
applied transformations transform G to the form given in equation (14). 
The encoded operators for this gauge group are clearly (Oifc|0) and (0|0/fc). 
We conclude with a simple example that illustrates the process. 

Example 17. To compare with the standard form method, we consider the 
[[4,1,1,2]] code again. Let the gauge group G, stabilizer S and encoded 
operators given by L. 
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In matrix form G can be written as 
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The transformations consisting ofT x = CNOT^CNOT^CNOT 1 - 4 followed 
byT 2 = I®H®H®H maps G to 
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Now transform the second row using T 3 = CNOT 2,3 CNOT 2 ' 4 . Then trans- 
form using T 4 = CNOT 4 ' 3 . We get 
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Applying T^ = H®H(^>I(^H gives us 
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We could have chosen T5 = if ® if ® I ® I , since i/te ejfeci 0/ if on £/te 
fourth qubit is trivial. The complete circuit is given as 
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Figure 9: Encoding [[4, 1, 1, 2]] code by conjugation method 

By switching the target and control qubits of the CNOT gates in T3 and 
T4 we can show that this circuit is equivalent to circuit shown in Figure 10. 
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Figure 10: Encoding [[4, 1, 1, 2]] code by conjugation method 

It is instructive to compare the circuit in Figure 10 with the one given 
earlier in Figure 4- The dotted lines show the additional circuitry. Since the 
gauge qubit can be initialized to any state, we can initialize \g) to j0), which 
then gives the following logical states for the code. 



|0) = |0000) + [1111) + |0011) + |1100) , 
|T> = |0000) + |1111) - |0011) - |1100). 



(15) 
(16) 



It will be observed that IIXX acts as the logical Z operator while IZIZ acts 
as the logical X operator. We could flip these logical operators by absorbing 
the H gate into \ip). If we additionally initialize \g) to |0) ; we will see that 
the two CNOT gates on the second qubit can be removed. The circuit then 
simplifies to the circuit shown in Figure 11. 
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H 



|0> 
|0) 



Figure 11: Encoding [[4, 1, 1, 2]] code by conjugation method - optimized 

This is precisely, the same circuit that we had arrived earlier in Figure 5 
using the standard form method. 

The preceding example provides additional evidence in the direction that 
it is better to initialize the gauge qubits to zero and avoid the encoding 
operators on them. 

Conclusions. In this paper, we have demonstrated that the subsystem 
codes can be encoded using the techniques used for stabilizer codes. In 
particular, we have considered two methods for encoding stabilizer codes - 
the standard form method and the conjugation method. While the standard 
form method explored here required us to initialize the gauge qubits to 
zero, it admits two two variants and seems to have the potential for lower 
complexity; the exact gains being determined by the actual codes under 
consideration. The conjugation method allows us to initialize the gauge 
qubits to any state. The disadvantage seems to be the increased complexity 
of encoding. It must be emphasized that the standard form method is 
equivalent to the conjugation method and it is certainly possible to use this 
method to encode subsystem codes so that the gauge qubits can be initialized 
to arbitrary states. However, it appears to be a little more cumbersome and 
for this reason we have not investigated this in this paper. There is yet 
another method for encoding stabilizer codes based on the teleportation 
due to Knill. We expect that gauge qubits can be exploited even in this 
method to reduce its complexity. It would be interesting to investigate fault 
tolerant encoding schemes for subsystem codes exploiting the gauge qubits. 

Appendix 

The logical states of a stabilizer code. We assume that our basis input 
states are of the form |0)® \a\ . . . a^), where aj G {0, 1}. Clearly, we have 
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freedom in the choice of the states into which each of these states are encoded 
to. Additionally, we have freedom in the choice of the encoded operators 
though they are not entirely unrelated. Perhaps, this is best illustrated 
through an example. Let us consider Shor's [[9, 1, 3]]2 code. A choice of the 
logical states for this code is 

|0) = (|000) + |111))(|000> + |111»(|000> + |111», 
|T) = (|000> — |111))(|000) — |111))(|000) — |m>). 

For this choice of the encoded states the logical Z operator is X® 9 and the 
logical X operator is Z® . On the other hand, let us see what happens if 
we choose the logical states as follows: 

|0> = 1000000000) + |000111111) + 1111000111) + 1111111000), 
|T) = llllllllll) + |111000000) + 1000111000) + 1000000111). 

In this case the encoded X operator is X® 9 and encoded Z operator is Z® 9 ; 
they are flipped with respect to the previous choice! 

So it becomes apparent that the assignment of the encoded operators as 
logical Z or X is flexible and it seems to depend on the choice of the logical 
states. But are we free to choose any basis of the codespace as the encoded 
logical states. We can show that this cannot be. For instance let us choose 
the logical zero state to be a superposition of the previous two assignments. 
Then we have 

|0) = (|000) + |111))(|000) + |111))(|000) + |111)) 
+ 1000000000} + 1000111111) + 1111000111} 
+ 111111000) . 

The possibilities for the logical Z operator 7 are ±X® 9 , ±Z® 9 , ±X® 9 Z® 9 . 
But for none of these operators we have Z |0) = |o). As these are the only 
possible encoded operators (modulo the stabilizer which acts trivially in any 
case), this is not a valid choice for |o). This raises the question what are all 
the possible valid choices for the logical states. Let us look at yet another 
choice of logical states. 

|0> = (|000> - |111»(|000> - |111>)(|000> - |in>), 
jl) = (|000) + |111))(|000) + |111))(|000) + |111)). 

7 Including scalar multiples of i will not change our conclusions. 
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In this case, the encoded Z and X operators are —X® and Z® respectively. 
This gives us a clue as to the possible logical all zero states for a given 
stabilizer code. The all zero logical state is the state in the code space that 
is fixed by the stabilizer and the logical Z operators. Assuming that S is the 
stabilizer and Cp n (S), its centralizer, we can can pick any k independent 
commuting generators in C-p n (S) \ SZ(V n ) as Z operators. Hence, we have 
the following lemma. 

Lemma 18. Let S be the stabilizer of an [[n,k,d]]2 stabilizer code. If L < 
C-p n (S) is any subgroup generated by n commuting generators such that L n 
Z{V n ) = I and S < L, then the state stabilized by L is a valid logical all 
zero state for the stabilizer code defined by S. 

The implicit choice of |o) made in Lemma 2 (by picking the encoded Z 
operators, at least the representatives) is convenient in the sense it allows us 
to speak of a canonical |o) without ambiguity. This |o) can be conveniently 
identified with the state P |0) , where it will be recalled that P is the 
projector for the stabilizer code given as 
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